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Abstract 

A systematic algorithm for obtaining recurrence relations for dimensionally 
regularized Feynman integrals w.r.t. the space-time dimension d is proposed. 
The relation between d and d — 2 dimensional integrals is given in terms of 
a differential operator for which an explicit formula can be obtained for each 
Feynman diagram. We show how the method works for one-, two- and three- 
loop integrals. The new recurrence relations w.r.t. d are complementary to 
the recurrence relations which derive from the method of integration by parts. 
We find that the problem of the irreducible numerators in Feynman integrals 
can be naturally solved in the framework of the proposed generalized recurrence 
relations. 
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1. Introduction 



Many phenomena in high energy, solid state and statistical physics can be described by 
quantum field theoretical models considered not only in four dimensional space-time 
but also in two-, three- and other space time dimensions [IJ. In many cases pertur- 
bation theory is the basic tool for calculating different physical quantities. The same 
Feynman diagrams occur in different models with different space-time dimension d. 
This parameter may be considered as a regularization parameter || different from the 
value of the physical space-time dimension. Usually one needs to set up a Laurent 
expansion of the diagram w.r.t. e — (21 — d)/2, with 21 being the dimension of the 
physical space -time in the problem under consideration. The coefficients of the ex- 
pansion are different for different I. The most advanced methods were developed for 
the evaluation of d = 4 — 2e dimensional integrals. For example, the standard package 
MINCER for calculating massless three-loop diagrams is now widely available and 
has been used in many multiloop calculations. 

Another reason to consider d different from its physical value was given in 
There, a new approach for obtaining nonperturbative information from a quantum 
field theory by expanding Green's functions as a series in powers of d was proposed. In 
light of this investigation one can try to formulate approximate methods for calculating 
individual Feynman diagrams at some asymptotic values of d, for example as \d\ — > oo 
or as \d\ — > 0. 

One of the most powerful methods for evaluating Feynman diagrams is the method 
of integration by parts || . In this approach dimensionally regularized Feynman inte- 
grals are considered as functions of the exponents of the scalar propagators. Integration 
by parts gives relations connecting integrals with some exponents changed by ±1, very 
similar to the relations for the contiguous hypergeometric functions. This is not a big 
surprise since it has been known for some time M that Feynman amplitudes belong to 
the class of hypergeometric functions. In many cases, integrals are just proportional 
to the well-studied hypergeometric functions 

Using the recurrence relations, a variety of Feynman integrals can be reduced to the 
restricted set of so called "master integrals" . The proof of completeness for the set of 
recurrence relations obtained by integration by parts and the problem of a systematic 
algorithm how to use them for an arbitrary diagram are open questions. Any additional 
information in this respect may be useful for solving both problems. 

In the present paper we propose a systematic formulation of the recurrence relations 
w.r.t. d. We shall show that these relations cannot be obtained by the method of 
integration by parts and therefore should be considered as an important addition to the 
Chetyrkin-Tkachov (CT) recurrence relations. As a concrete example, we demonstrate 
that CT recurrence relations with the recurrence relation w.r.t. d compose all possible 
recurrence relations for the considered integral. New recurrence relations are not of 
purely academic interest. We demonstrate how one can calculate Feynman integrals 
using these relations. We also rederive some useful relations for one-loop diagrams with 
arbitrary number of external legs. The relations connecting integrals with different d 
may be also useful for calculating integrals with d ^ 4. 

We show how tensor integrals can be represented in terms of integrals with the 
changed d. This representation allows us to write Feynman integrals with irreducible 



numerators as a combination of scalar integrals having different values of d. Thus, the 
solution of the generalized system of recurrence relations automatically leads to the 
solution of the problem of irreducible numerators. 

It turns out that the new recurrence relations are especially useful in the new 
method for the momentum expansion of the scalar Feynman integrals proposed in [ |T0|j . 

The paper is organized as follows. In Sect. 2, we present the main ideas of our 
method. First, using the parametric representation we derive relation for arbitrary 
scalar integrals with different d. Then we show how tensor integrals can be expressed 
in terms of combinations of scalar integrals having different values of the space-time 
dimension. In the framework of this approach the solution to the problem of irre- 
ducible numerators in Feynman integrals is proposed. An algorithm for obtaining new 
generalized recurrence relations including integrals with different d as well as different 
exponents of the propagators is developed. 

In Sect. 3, explicit relations connecting integrals with shifted d for some one-, two- 
and three-loop integrals are derived. 

In Sect. 4, analogous relations for the one-loop integrals with an arbitrary number 
of external legs and arbitrary powers of propagators are given. We reproduced already 
known results for these integrals and obtained more general new ones. 

In Sect. 5, we demonstrate how one can calculate integrals explicitly by solving 
recurrence relations w.r.t. the dimension of space-time. 

In Sect. 6, the new recurrence relations are used to calculate the e = (6 — d)/2 
expansion for the two-loop self-energy diagrams from the expansion developed at e = 
(4-d)/2. 



2. Relations for integrals with different d 

The subject of our consideration will be dimensionally regularized scalar Feynman 
integrals. An arbitrary scalar L loop integral can be written as 

L N 

G {d) ({si},{™ 2 s }) = IlJ d%YlPg mj , (1) 



i=i 
where 

-y L E 

Pk,m = 7~j~a ™2 i = Y U jnkn + Y VjinQm, (2) 

[K - m -+- 1€) n=1 m=1 

q m are external momenta, {s«} is a set of scalar invariants formed from q m , N is the 
number of lines, E is the number of external legs, to and r\ are matrices of incidences 
of the diagram with the matrix elements being ±1 or (see, for example, Ref. [|TT|j). 

To find the desired relation we shall use the parametric representation of the inte- 
gral, which can be found in the literature [[11]], Jl2|]. For an arbitrary scalar Feynman 
integral in d dimensional space-time we have: 

G«({«>.{m:» =f (-fn ^ r ■ f^C e«-£^.«-,>, (3 ) 



where D(a) and Q({si}, a) are homogeneous polynomials in a of degree L and L + 1, 
respectively. They can be represented as sums over trees and two-trees of the graph 
(see, for example, fHj): 



d{oc) = e ( n -",-)• 

over over 
trees chords 

Q({si},a) E II ( E <?) 2 - (4) 

of er of er over comp. 

2 — trees chords of 2 — tree 

These polynomials are characteristic functions of the topology of the diagram and of 
its subgraphs. Since D and Q will play an important role in the rest of the present 
paper we remind the reader of the definitions of the trees and two-trees for connected 
diagrams. Any connected subdiagram of the diagram G containing all the vertices of 
G but is free of cycles (loops) is called a tree of G. Similarly, a two-tree is defined as 
any subdiagram of G containing all the vertices of the original diagram, but is free of 
cycles, and consisting of exactly two connected components. Finally, a chord of a tree 
(two-tree) is defined as any line not belonging to this tree (two-tree). Examples of D 
and Q for some particular diagrams will be given in the next sections. 

In the case when the z/,-'s do not depend on d, one can see from (^), that d enters 
the integrand in a rather simple way. Only the exponent of D(a) linearly depends on 
the dimension of space-time. 

In order to find the relation between integrals in different dimensions d we have 
to assume at first, that all scalar propagators in (^) have different masses. Next, we 
introduce the polynomial differential operator 

which is obtained from D(a) by substituting ctj — > dj = d/dra 2 -. The application of 
the operator D(d) to the integral (|3|) gives D(a) in the numerator of the integrand: 

D(d)e- i ^ airn ' -f J D(a)(-z) L e"^ aim '. (6) 

The resulting integral is proportional to the same integral with d changed to d — 2: 

G( d - 2 )({s,},K}) = (- i)^ D(d)G^({ Sj },{m 2 s }). (7) 

After having performed the differentiation we identify masses with the one's of the 
original integral. 



We may include tensor integrals into our consideration using standard methods [JL1 
13]. To each line one introduces an auxiliary vector aj and then differentiates w.r.t. 



these vectors. The parametric representation for a tensor integral with products of n,- 



vectors corresponding to the j-ih line reads: 



L . N ni n N , , i>± N A-Vj-rij 

n / 4% n % n k, n = * L (f nW 

i=l J j=l 1=1 s=l Vi/ j=l L \ U ]> 

r =i da lflr ^ da NXa Jo Jo [D(a)]? 



where 



™f = ™i + ^2, (9) 



Sj are scalar invariants formed from vectors q i 



1i = ( li + Y, e ij a i^—^ ( 10 ) 

the qiS are the external momenta incoming at a vertex i and e is the incidence matrix 
defined as: 



+1 
-1 




if the oriented line j points away from the vertex i 
if the oriented line j points toward the vertex i 
if the line j does not contain the vertex i. 



Differentiation w.r.t. aj will produce external momenta and metric tensors times 
some polynomials R s (a) divided by D(a) to some power in the integrand. The poly- 
nomials R s (a) have to be converted into operators R s (d) and the powers of D(a) are 
absorbed into the redefinition of d. In this way any tensor integral will be expressed as 
a sum over a set of tensors formed from external vectors and metric tensors multiplied 
by a combination of scalar integrals with the shifted value of d. At the one-loop level 
such a representation was already proposed in IfUfl . 

Tensor integrals in momentum space can be written in terms of scalar ones without 
direct appeal to the parametric representation (§). The procedure described above 
may be derived from the following formula: 



L „ N ni n N L . N 

n / d% n % nj n ^ n w = n q , o, d + } n / d% n ^ w 

8=1 j=l r=l s=l i=l j=l 

where the tensor operator T is: 

p -iQ{{si},a)p 

n q ,d,*+)= . ni+ ^ 

„2 



x fi — n -^—e i[Q{{ ^ a) -^ ^ o{a)]p 



Qj = 

otj = idj 



;i2) 



Here d + is the operator shifting the value of the space-time dimension of the integral 
by two units: d + G^ = G^ +2 ^. As before we have to assume that at the beginning 



all propagators have different masses and after applying the operator T to the integral 
one should set masses equal to the required ones. 

The representation ( |TTD for tensor integrals may give a solution of the problem of 
irreducible numerators, i.e. the appearance of scalar invariants in the numerator which 
are absent in the scalar propagators. This kind of integrals can be expressed in terms 
of scalar integrals without numerators but with the changed space-time dimension. 
Therefore the combined set of recurrence relations, i.e. relations with the changed d and 
those obtained from integration by parts, should be used to reduce all integrals to the 



set of scalar master integrals. In |K| we presented a method for the small momentum 
expansion of multiloop scalar integrals based on (^-recurrences. This method from 
the very beginning does not produce integrals with irreducible numerators although 
those appear in traditional methods for the small momentum expansion. In principle 
irreducible numerators can be regarded as propagators raised to negative powers and 
the corresponding integrals can be considered as object in a more complicated class of 
integrals with additional new denominators. In our approach one remains in the same 
class of functions, satisfying generalized recurrence relations i.e. recurrence relations 
derived from the method of integration by parts plus relations connecting integrals 
with different d. 

Having representation (|TTD at hand, we can now state the procedure for obtaining 
new, generalized recurrence relations. The starting identity: 



tj**£z fe^fi'SJ-O. (13) 



d_ 

where Xi are arbitrary constants, written in this form turns out to be rather convenient 
for the derivation. After performing the differentiation one would usually express 
scalar products with integration momenta in terms of invariants that occur in the 
denominator of the propagators. At this point we propose to change the derivation 
and write all or only some of integrals containing scalar products with loop momenta 
as a combination of integrals with changed d. In this way we produce many relations 
including integrals with changed exponents of scalar propagators and changed values 
of the space-time dimension. Combining the different relations (or choosing in a proper 
way xi in ( pf) ) one can try to find the most optimal set of relations for the reduction 
of the concrete class of integrals to the set of basic "master" integrals. In fact, the 
method of integration by parts corresponds to some specific representation of scalar 
products in (P~3"|). Our derivation is more general and it includes integration by parts 
method as a particular case. 

We expect that solutions of the generalized recurrence relations and those obtained 
by the method of integration by parts will be of the same complexity. The fact that in 
case we know the explicit result for an integral in terms of hypergeometric functions, 
d/2 appears in the same manner as any exponent of the propagators, can be considered 
as a confirmation of this statement. However further investigation of this problem is 
needed. 

Several remarks are in order. All relations which connect integrals with shifted d 
are valid for arbitrary momenta and masses and also for their real and imaginary parts. 
We want to stress here that in general scalar products in (|7|) must be considered as 



independent variables. One cannot use any restrictions valid for some specific values 
of d. 

All relations obtained with the help of the parametric representation may be prof- 
itably used in the frame of momentum as well as configuration space. 

The minimal change of the space-time parameter is two. From concrete examples 
0"0) when the result is known in terms of hypergeometric functions, one can observe 
that changing d by ±2 we obtain contiguous functions and therefore may hope to 
find relations between integrals. The change of d by ±1 shifts the parameters of the 
hypergeometric functions by ±1/2 producing functions which does not belong to the 
class of contiguous functions. In general there are no relations between those functions. 

Application of the differential operators to the integral will increase the powers of 
the denominators. To simplify the r.h.s. of (0) the CT recurrence relations can be 
used. In general every integral in d — 2 dimensions can be reduced to combinations of 
a rational function of scalar products of external momenta, masses and d and rational 
multiples of master integrals 2^ ({sj}, {m 2 s }): 

^(W.K^EQfel.Kl.^f ({*},{*#)■ (i4) 

The set of basic integrals in d — 2 dimensions can be expressed in terms of the same 
integrals in d dimensions: 

^ a) (W,W}) = E s «(W.W},^i W (W,{m2}). (is) 

3 

This relation can be inverted and therefore we will get a representation of an arbitrary 
d dimensional basic scalar integrals in terms of d — 2 dimensional ones. Such a rep- 
resentation may be useful for some practical calculations. The most evident practical 
application of the relation connecting integrals with different d: having the e expansion 
for some particular do we can find similar expansion in do ± 21 (/ integer) dimensions. 

The proposed relations can be used, for example, in the evaluation of massless 
propagator type 2 — 2e dimensional integrals using the package MINCER || writ- 
ten for 4 — 2e dimensional integrals. The polynomials Q({si}, a), D(a) can be easily 
constructed by means of a computer program for any particular integral |i"5|| . 

3. Examples 

In this section several illustrative examples of the new recurrence relations will be 
presented. We start with the one-loop propagator type diagram with massive particles: 

« 2 (? 2 , ml, ml) = J 0^P^ mi P^ q ,mr (16) 

In this case D(a) = a± + a 2 and therefore 

ll d - 2 2) (q 2 ,mlm 2 2 ) = -(d, + d 2 )4%{q 2 , mf, m 2 2 ) 

= U2 (q 2 ,ml,ml) - v 2 I^f „ 2+1 (q 2 , ml, ml). (17) 



Here and in the following a,i is assigned to the line (or propagator) with mass m«. We 
can get another recurrence relation connecting integrals with different d following the 
prescription given in Sect. 2. From the identity: 

/ + q ^ P ^n P kU m2 ] = 0, (18) 

we obtain (omitting for the moment arguments of /^L, 



d d ki 1 

Vl J J~d/^( qk ^ P kl mi P ki-q,Tn 2 

= (| " ^ 4% - ^4f-i , 2+ i - V2 + »M - ml)llf V2+1 . (19) 

The integral with the scalar product (qki) can be written as a scalar integral with 
shifted d. According to (Q) and ( |TUD the function Q with auxiliary vectors ai, a 2 , 
which are related to the lines with masses m 1; m 2 , respectively, is: 



Q({s},a) =[q + 



ai a 2 



2ai 2a 2 

Formulae (|TT]) and (0) for this case yield the relation: 

d d h 



a\a 2 . 



(20) 



r d fci +1 _ 2 ( rf+2 ) 

7 [iTfd/2] ^ ^ r fe,mi r «!i- g ,m 2 ~ ^2<? 1/3+1 • 



(21) 



Inserting (|2l|) into (0) we obtain the desired identity: 



WiffiUi " (5 "^i)/S 2 + ^ „ + i + ^i4?+i V2 -Hf- m?)4?„ 2 +i= 0. 

(22) 

In addition to the above relations two more relations can be obtained from the tradi- 
tional method of integration by parts: 



2v 2 m 2 l i 



2 Ad) 



V\ U2 + 1 



+ fi/Jf-i „ 2 +i + ^iK + m 2 2 - q 2 )ll d) +1 V2 -(d- 2v 2 - ujW = 



^4i+i ^-i-^^-i i , 2 +i-^i(^?-^2+9 2 )4 1 +i ^-^(rJ^-ml-g 2 )/^^! 

+ (^ 2 -^)/£ 2 = o. 



(23) 



For the simplest case, at mi = m 2 = and V\ = v 2 = 1 from (pT7|) , (|23"D we readily get: 



/{TV, 0,0) = -2/£V,0,0) 



2(d-3) 



A ( iV,o,o). 



(24) 



This formula can be easily verified from the explicit result for lj)fl 2 (q 2 , 0, 0). 



At mi = 0, m 2 = m, the relations connecting integrals lj)fl 2 (q 2 , 0, m 2 ) with different 
d are: 

W/^+iC? 2 , 0, m 2 ) - (± - J /£ 2 (g 2 , 0, m 2 ) + x , 2+1 (g 2 , 0, m 2 ) 



-^g 2 if, 2+ i(g 2 ,o,m 2 ) = o, 

ll d ~ 2 2) (q 2 ,0,m 2 ) + , 2 (g 2 ,0,m 2 ) + v 2 I^ ^(q 2 , 0, m 2 ) ee 0. (25) 

The integral 7^ 2 (g 2 ,0, m 2 ) is proportional to the Gauss hypergeometric function 0: 



(a 2 



m 



i>i + 1/ 2 - |)r( 



"1, 



(m 2 y i+u *-2T(l)T(is 2 



2 r l 



vi, vi + "2 - f ; <? 



d . 

2 ' 



(26) 



As is well known there are fifteen relations of Gauss between contiguous functions 2 Fi • 
Substituting (p6|) into (p3|) one can find correspondence between the CT recurrence 
relations and only six relations of Gauss. The reason is obvious - in the CT relations the 
third parameter of 2 Fi in ( ^6|) does not change and therefore all corresponding relations 
for contiguous functions cannot be reproduced. If we include into consideration also 
identities (p5|) we cover all fifteen relations, though in principle, to reduce the integrals 
J^ 2 (g 2 , 0, m 2 ), with integer v\ and v 2 to two boundary integrals, it is enough to know 
only the CT relations. 

Now we consider two-loop bubble integrals with three different masses: 



The function D(a) for this integral is: 

D(a) = (X\(X2 + ai«3 + a 2 a^, 

and hence, 

Jnlr = (^12 + ^13 + ^23 ) ^V X V2.V 3 ' 

In the above formula and in what follows we use the shorthand notation: 



(27) 

(28) 
(29) 



0\ 



ii...i N 



N Q 

TT — — 



Let's take for simplicity V\ — v 2 — v% — 1. In this particular case the relation connecting 
J111 with different d first was found in [16|]. Exploiting the relation obtained by the 
method of integration by parts: 



A diJ { n\ = J(d - 3)(fcA) Jffi + (d - 2) J&{ + 



r(«0 



(d-2) r 
4m 2 



(ft A) Jg} + (ftA)J^ , (30) 



where 



a /I /I /l OO 00 o 

A = m 1 + m 2 + m 3 — 2m l m 2 — 2m 1 m 3 — 2m 2 m 



3- 



(31) 



and similar relations for c^Jm, ^3 J11L from fl29|) we reproduce the relationship between 
integrals with different d given in |L6| : 



(d-2) 



11 



( d - 2)(d - 3) jffl + ^-T 2 I 2 - - ) [ 8l82 d 3 A + Sl s 3 d 2 A + s 2 s 3 d l A], (32) 



d 



where s$ = 

One can easily obtain another recurrence relation for the integral (|27|) with different 
d. From the identity 



d d h d d k 2 



d d 
+ 



s dk^ dk 2il/ 
keeping the scalar product (kxk 2 ) untouched, we get: 



1^ fci,mi k\—k2,m,2 k2,rrt3 



= 0, 



(33) 



71° 



W2 

k-i—k2,m2 ± k2,m% 



P"3+l i j. m 2 Ad) 



(34) 



The first integral in (0) can be expressed in terms of integrals with another value of 

ion d by i 

Q({s},a) 



the space-time dimension d by using formula (12) with 



a 3 \ 2 
4 \oi ct 2 a 3 J 



a 2 
ol 2 



(35) 



Here again vectors dj correspond to lines with mass mj. In order to obtain the integral 
with the scalar product (kik 2 ) in the integrand one has to differentiate w.r.t. a\ and 
a 3 which leads to the relation: 

1 f 1 d 2 

— y d d A;i d d k 2 {k x k 2 ) P^ <mx P^-^ m2 P^ mz = -u 2 dJ Vl „ 2+1 „ 3+1 . (36) 
Substituting fl3"6| ) into (|31) we obtain: 

Wdj^+i ,3+i + 2^i JS-i , 2 , 3 - - 2^i)41,3 = 0- (37) 

This identity was used to evaluate the coefficients in the small momentum expansion 
of the two- loop master diagram with all masses equal by the method proposed in [[nj . 

At the three-loop level, bubble diagrams in general have the topology shown in 
Fig.l, where each line corresponds to a scalar propagator with an arbitrary exponent. 




Fig.l. Three — loop bubble diagram 

The relation between d and d — 2 dimensional three-loop vacuum integrals with all 
masses arbitrary is: 



G {d - 2 \{m 2 s }) = -i + d\ 2A + d 3 126 + df 35 + «9 3 36 + d 3 U5 + + & 

7T J v 



3 

156 



+d 3 2U + 4* + <9 2 3 45 + dl m + «9 3 56 + «9f 45 + dl m + «9 3 56 ) G^({m 2 s }). (3* 



The numbering of the masses corresponds to the numbering of the lines in Fig.l. 

Three-loop master bubble integrals encountered in the small momentum expansion 
of the QED photon propagator were studied in . Every integral in this case can be 
expressed in terms of three basic structures. The only nontrivial combination of two 
integrals taken as one of such structure was: 



m 3d - 12 B[ d) 



1 



(d-2)(d-3) 



d d h d d k 2 d d k 3 



I)] 3 



P P P 

1 k2,m 1 k-i-pi 1 k2— ^3,0 



x [-ffci.m -Pfci-ifc3,0 Pk 1 -k 2 ,0 ~ Pki,0 P) 



k\,0 -Oci— fe3,m Pk\—k2,m 



(39) 



By using (|38|) and the method of integration by parts we find the following relation for 



B 



(d). 



B 



(d- 



_ 2) _ 3(d - 4) 2 (d - 5) (3d - 14) (3d - 16) {d) 33d 3 - 367d 2 + 1170d - 864 

Ba -\- 



16(d - 6Y 



8(d-4) 2 (d-6) 5 



7(8d 2 - 80d + 195) T(f - l)r 2 (6 - d)r(9 - f ) 

(d-4) 2 (d-6) 5 r 2 (3- f)r(i2-2d) ■ 



(40) 



We observe that the relation includes B\ itself and two additional terms with a trivial 
d dependence. The results for other diagrams of the three-loop photon propagator look 
similar, i.e. they exhibit three terms of a structure like in (|40|) . 

Other useful three-loop vacuum integrals were introduced in [IT8J: 



m M - 12 D 3 



rn 



3d-12 



B, 



d d k x d d k 2 d d k 3 

[m d / 2 r(3 - d/2)] 3 

d d k 1 d d k 2 d d k 3 



d / 2 r(3 - d/2)] 



\17C 



Pk\,oPk2fiPks,mPk\—kifiPk\—k2,mPk2—kz,mi (41) 

(42) 



_ p p p p 

g 1 kifTn 1 k'i^m 1 ki—k2,m 1 k2— ^3,0 



These integrals occurred in the evaluation of the three-loop QCD correction to the 
electroweak p parameter. The space-time recurrence relation for looks more com- 
plicated than that for S4 and it includes also B^: 

(d _ 2 ) _ 16(63d 3 - 832d 2 + 3622d-5176) 64T(| - 1)T(5 - d) 



D 



(d - 4) 4 (d - 5) 2 (d - 6) 3 (d - 4) 3 (d - 5)(d - 6) 3 T(3 - £ 



x 



1 + 



(7d-32) r(f -l)r(7 
3T 2 (3-f) 



+ 



(37d 2 - 350d + 828) T(5 - d)r(7 - 



2 ) 



24(2d-9) r(3 - f )T(9 - 2d) 



4(d - 2)(d - 3)(3d - 8)(3d - 10) 
9(d-4)(d-6) 3 



4(d-2)(d-3)(d-4) 



(d-6) ; 



D 



(43) 



The corresponding relation for £5 is simpler and reads: 



B 



(d-2) _ _4(d - 2)(d - 3) (3d - 8) (3d - 10) 
5 ~ 9(d-4)(d-6) 3 



B 



(d) 



64(15d- 52) 
9(d-4) 4 (d-6) 3 ' 



(44) 



It can be obtained from (|38|) by observing that c^G^ = d^G^ = 0; which holds 
because corresponds to the diagram Fig.l with shrinked second and fourth lines. 



For we found an explicit expression in terms of hypergeometric functions 3F2 and 



2F1 with the argument 1/4, satisfying relation ([441). 

Relations (^0D,(^3|) and ( f44|) can be used for the computation of the coefficients in 
the small momentum expansion of Feynman diagrams by the method proposed in [|Kj . 



A more detailed analysis of the generalized recurrence relations for two- and three- 
loop diagrams will be given in a future publication. Several examples of these recurrence 
relations one can find also in 



4. Relations for n-point one-loop integrals 

In this section we consider scalar one-loop integrals depending on n — 1 external mo- 
menta: 

where 

bl = (<? — Pj) 2 — ntf + it for j < n and [n] = q 2 — + ie. (46) 
The corresponding diagram and the convention for the momenta are given in Fig.2. 



q-pi 




Fig.2. One-loop diagram with n-legs 
For the diagram Fig.2 the function D(a) is a polynomial linear in af 

n 

£>(«) = X>* (47) 

and hence 

n 

iL d - 2) = -Y,dA d) - (48) 

3=1 

To get rid of the derivatives in this relation we use the method of integration by parts. 
Let's consider the identity: 



that is valid for arbitrary Xj. Upon differentiation and expressing scalar products in 
the numerator in terms of invariants [j] in the denominator, we get: 



f ill n \ n 1 / n_1 v 



1/,- n 



= 0, 



(50) 



where R) n \{xi}) = ££=i with 

2 2,2 

Rif = mj - m l - p\ + 2 PiPj, 



R: 



(») 



m i — m n ~ Pi f° r i < n , 

for i,j<n 



(51) 



The integral related to the last sum in fl5U| ) can be made proportional to the r.h.s. of 
(f48D - This requires to find Xj such that 



= 1, j = l,...,n. 



(52) 



Imposing these conditions we obtain a system of n equations for Xj. The solution of 
the system reads 



Xi = - di In A n for i < n, 



where G n is the Gram determinant 



Gn-1 

2A n ' 



(53) 



VlPl P\P2 ■■■ V\Pn 
V2V\ P2V2 ■■■ V2Vn 



VnPl VnV2 ■■■ PnPn 

and A n is proportional to the determinant 



of the n x n matrix C: 



-2^et(C), 



(54) 



(55) 



where it is assumed that m = m n and po = 0. Substituting Xj into (|50|) , we obtain 
the following relation between d dimensional and (d — 2) dimensional integrals: 



n-1 



- v n ) (djA n ) + (d-u n -J2 Vj)G n -i 

j=± 0=1 



#° = E id An) InT ] + 2A ^~ 



■2) 



(56) 



The index j of 7^- means that the factor in the integrand must be changed 

into [j]/^]^- For Uj = 1 we obtain the simpler relation: 

n 

( d - n - l)G n _ 1 7W = £ (d,A n ) 7^ 2) + 2A„7( d " 2 ). (57) 

i=i 

For d = 6, n > 6 and assuming that external momenta are four dimensional the 
Gram determinant G„_i vanishes and hence the term with 1^ drops out and we get a 
relationship between 4-dimensional integrals with n legs and integrals with n — 1 legs. 



Such a relation was first obtained in (see also |20|). When d = 6 and n = 5 the 



l.h.s. of ( |5TD is zero and we arrive at a formula for reducing 4 - dimensional pentagon 



integrals to box integrals. This relation was first derived in [21]. For arbitrary d and 
n the equation (|57|) was obtained by a different method in Eq. fl56|) as far as we 
know is new. 

It is evident that by using fl7|) and integration by parts one can also derive relations 
similar to ([57]) for multiloop integrals. 

Equation fl56|) can be used in the reduction of the one-loop tensor integrals to scalar 
ones. In (JTJJ an explicit general formula for one-loop tensor integrals was derived: 



n- =1 [j? 



Mi— Mm 



A, fci, . . . , 1 

2A + E^i = M 



. M+2M-2X n-1 ( v .\ 



where (v)k = T(z/ + £;)/T(z/) is the Pochhammer symbol. The shorthand notation 
{[g] A [pi] fcl . . . [p n -i] fcn_1 }^i...^ Af corresponds to the symmetrical (w.r.t. \i\ . . . Hm) tensor 
combination, each term of which is constructed from A metric tensors g, ki momenta 
Pi, . . ., k n -i momenta p n -i- For example, 

In fl5"8|) A, ki > 0, max fc; = M, max A = [M/2] (integer part of M/2). For more details 
see ||14|| . Thus, the procedure of calculating one-loop diagrams will be as follows: one 
should use first (|58D , perform contractions of indices and then using ( |56|) to reduce 
all scalar integrals with the changed dimensions to the d = 4 — 2e dimensional set of 
integrals. After that scalar integrals should be reduced to the set of master integrals 
by making use of recurrence relations obtained by the method of integration by parts. 



5. Evaluation of one- loop integrals using space-time 
recurrence relations 

With the help of equations (|TJ|) and one can obtain relations that include only 
one particular integral with different shifts in d: 

E **({*,}, d)I^ k \{ Si }, {ml}) = I ({s t }, K}, d), (59) 

k 



where Jo({sj}, {m 2 }, d) is some explicitly known expression. It can be eliminated from 
the equation giving rise to higher order recurrence relations. For an explicit solution 
equation (^) looks simpler than the ones obtained by the method of integration by 
parts. In this section we will demonstrate how the space-time recurrence relations 
proposed in the present paper can be used for evaluating integrals. 
We consider the following one-loop propagator type integral: 

jM = m 4-^ f J^L_ P ki>m P (60) 



By using (0) and the first relation of fl23|) we obtain: 

,2 orfo _ d\^2 



/(d _ 2) = 2(d-3)m 2 l(d) _ 2T(2 - f )m 2 
(q 2 — Am 2 ) (q 2 — 4m 2 ) 



(61) 



Introducing 



1 \4m 2 -q 2 ) r(f-l) ' [b ] 



we get the simpler equation: 



m 



r(d-3)r(2-|) 



F(d-2) _ F(d) _ ^Hf 2^ ffi oN 

^4m 2 -?v ^r(f-i) ■ 1 j 

Without loss of generality, one can parametrise d as d — 21 — 2e, where I is an integer 
number and e is an arbitrary small number. The inhomogeneous term in ( p3f) can be 
absorbed by a redefinition of 7^: 

2 £ - 3) r(2- J + £ ) +c , 




i2i-2 £ p(j _ 1 _ e ) 

yielding the following very simple equation for C\: 

C[ = C 1 - 1 . (65) 

Since / in our case is an arbitrary integer we can conclude that C[ does not depend 
on I at all. It can be found, for example, by taking the limit I — > oo, or |d| — > oo. 
Taking the limit \d\ — > oo is quite a delicate matter. The larger \d\, the more divergent 
becomes the integral. In order to keep the regularization of the integral, we have to 
consider e as complex and I large. The asymptotic value of the integral as I — > oo can 
be found by the method of steepest descent. It depends on the sign of q 2 . From the 
parametric representation of the integral one gets for q 2 < and / — > oo : 



—Tim \ / g x 



/( .0^ r(2 _ i + £) (_^L) (1 + 0(1/0). (66) 

At large Z and g 2 < 0, the sum in ( p4|) is convergent which allows us to find C e : 

i 2£ n m 4 ~ / m 2 V" £ r(2j - 2e - 3)T(2 - j + e) 



C" 

£ sin vt£ v ^ 2 "(4m 2 - g 2 ) 3 / 2 ^ V 4 ™ 2 ~ ?V i 2j - 2£ T(j - e - 1) 



(67) 



Substitution of C £ in the above formulae yields: 



Vir(j-l)r(2 
r (¥) 




d-i . 

2 ' 
d . 

2 » 



4m 2 



4m 2 — g 2 



(68) 



By using the formula for the analytic continuation of the hypergeometric series we get 
the known result: 



/(<*) = r 2 



1,2- 

3 . 
2 ' 



2 > 



Am 2 



(69) 



The example considered illustrates the main ideas of how to use the relation (|7D for 
the evaluation of Feynman integrals. The same method applies without modification 
to more complicated cases. For example, to the one-loop integral (|T6|) with different 
masses at V\ — v-i — 1 or to the two-loop integral Q2"7| ) with nu\ = z/ 2 = u 3 = 1. The 
recurrence relations in these cases are very similar to (|6T|). We applied this technique 
also to some two- and three loop integrals. The main difficulty encountered in these 
computations was to find the asymptotic value for the integral for large d. An expansion 
at large d is frequently used in solid state physics and statistical physics. Our experience 
shows that it can also be used for the approximate evaluation of individual Feynman 
integrals. Details of these calculations will be given in a future publication. 



6. 



The two- loop self-energy for the 3 model in 6 — 2e 
dimensions 



The aim of this section is to show how <i-recurrence relations can be used to obtain the 
e expansion of the nontrivial integrals in one dimension from the known expansion in 
another dimension. 

We consider two-loop self-energy diagrams which are encountered, for example, in 
b 3 model with two massive and one massless scalar fields 1123 1 



the 



At the two loop 



level there are two nontrivial diagrams, given in Fig.3, contributing to the self-energy 
of the massless field. 





Fig.3. Two-loop scalar diagrams with massive loops 

The solid lines in Fig.3 are related to massive particles and the dashed lines to massless 
ones. The corresponding integrals read: 



n 



(d) 



d d k 1 d d k 2 



Z7rlr(4- f) 



P P P 



ki ,m^'k2 ,mf~k\. — q,m,P. k2 — q,m,I~k\ — ki ,0 ? 



i.Pk 



Using the e expansion in d = A — 2e dimensions we will find the e expansion of 
these integrals up to 0(e) at d = 6 — 2e. This can be done in three steps. 

The first step consists of expressing the integrals l[ d ' and in terms of basic ones, 
chosen in || to obtain the d — *■ 4 limit for the two-loop photon propagator. In fact two 
integrals from this basis are the integral l[ itself and its derivative w.r.t. the mass. 
The remaining basic integrals were the derivative w.r.t. the mass of the one-loop scalar 
integral with two massive lines, its square and a one-loop vacuum integral with a mass. 

Thus, at this stage we need to compute only the integral Ii- By using the recurrence 
relations given in || for integrals contributing to the two-loop photon propagator, we 
get: 



r(d) 



m 



(d-3) 



(d-3) 



1 - 



Am 2 



V 1 ) + (d-? 



2d- 10 



(d-3)(d-4) 



(71) 



where the prime denotes differentiation w.r.t. m 2 and 77^ is a basic integral occurring 
in one-loop calculations: 



(d) 



2 7\ 



1 Q 2 • 

1, O 2 , 

3 . 
2 ' 



q 



Am 2 



1 



(d-6) 
d-5 ' (d-5) 



+ 



1 



Am 2 



(d-2) 



(72) 



The next step is to find relations between the basic integrals l[ , 
d — 2 dimensions. Using Eq.(|7|), with: 

-D(a) = OL*,{OL\ + «2 + «3 + CK4) + («1 + «3)(«2 + CK4) 



7^) in d and 



(73) 



(en's are labeled as shown in Fig. 3) we get a system of equations in the form of fll5|) . 
This system enables us to find the required relations: 



3z{d - 3)(d - 4) 2 (3d - 14) (3d - 16)if } = — (d - 8) 2 (d - 6) -(d - 6)(d - 4)* 



+12(d - 4) (2d - 11)^ + 12(d - 7) (3d - U)z - 16(d - 7) (2d - 9) 7} 



(d-2) 



m 



+ — (d - 8) 2 (d - 6) [2(d - A)z s + 2(3d - 16)z 2 - 8(9d - Al)z 



+32(2d-9)] (i^)' 



2m 



2d-10 



(d-6) 2 (4-^ 
8(d - 4)(19d 2 - 188d + 462)^ 2 + 16(22d 3 



(d - 4) (23d 2 - 226d + 552)^ 3 



+64(2d-9)(d-5) 2 ] (H[ d) )' 



+ 



16m 



2c2-10 



(d-6) 2 (4-^) 



310d 2 + 1447d- 2238)^ 



(19d 2 - 182d + 432)^ 2 



-(56d 2 - 552d + 1344)^ + 16(2d - 9)(2d - 11) 77} 
16m 2d-io 



(d) 



(d-6) 2 (d-4)(d-5)(4-^) 



(d - 6)(d - Afz 2 - 2(d 3 + 3d 2 - 98d + 2i 
-8(d-5)(5d-28)(2d-9)], 



(74) 



-3z(d - 3)(d - 4) 2 (3d - 14)(3d - 16) (if )' = ^(d - 4)(d - 6)(d - 8) 2 

2 



X 



(d - 4)(d - 5)z 2 + (7d - 34) (d - 7)z - 4(2d - 9)(d - 7) JJ 



(d-2) 



+ ^-(d - 4)(d - 6)(d - 8) 2 [(5d - 24> 2 - 4(7d - 33)z + 16(2d - 9)] (if ~ 2) ) 



(d-2)V 



+ 



24(d - 4)m 2d - 12 
(4-z)(d-6) 2 



(d - 4)(7d 2 -70d + 174)^ 2 - (30d 3 - 424d 2 + 1986d - 3084)z 



-8(2d-9)(d-5)'](^-^ 



2 96(d-4)m M - 12 



(4-z)(d-6) 2 



[(d-4)(d-6)z 



4f2d 9V2d 11)] I 48(7d-32)(d-6>-192(2d-9)(5d-28) 
-4(2d-9)(2d-ll)j// 1 + m ^{A-z){d-Q) 2 ' (75) 

where z = q 2 /m 2 . 

The final step is to perform the e expansion up to a constant term. It is remarkable 
that the most complicated integrals l[ d ~ 2 ^ and (j[ d ~ 2 ^ at d = 6 — 2e give contribution 
starting; from the first order in e and therefore may be disregarded in the considered 



approximation. To calculate and I^' up to a constant term in e we need to expand 
to second order. The required expansion is: 



rid) 



(6-2e) 



1 + e 



2 + -ln(u) 



1 - 2/ x 2 , , . l l2 

- lir 



4 + — + — ln 2 (t;) + -hnV) 
6a 2a a 



1 + a 



Li 2 



,(76) 



where 



q 2 — 4m 2 ' 
from which we obtain 



v = 



1-a 
1 + a 



and Li 



i 2 (x) = - / 



* ln (l - y)dy 



4e _ 2 r(6-2e) 



3 -a 2 



+ 



31 - 7a 2 



(77) 



19 

+ 77 ~ 



+ 



47 



2h(a) 



6£ 2 (l-a 2 ) 12e(l-a 2 ) 24 12a 2 6(1 - a 2 ) 3a(l - a 2 ) 



5 



1 



7 1 

+ 



1 



+ 



48 3a 24a 2 48a 4 3(1 - a) 
48a 



ln» 



43 1 

- -7 + 



12 V a a 3 (1-a 2 ) 



\n(v), 



(78) 



m^- 2 lt 2£) 



15 + 28a 2 



+ 



389 



1207 



18e 2 (l-a 2 ) 108e(l-a 2 ) 324 36a 2 648(1 - a 2 ) 



l-3a 2 , . , If 4 

— — 7- rr a(a) + 777 10+- 

18a(l - a 2 ) v 7 144 \ a 



(1-a), 



\n 2 (v) 



+ 



43a 



^54a 36a 3 54(1 -a 2 ), 



ln(v), 



(79) 



with 

h(a) = In 2 (1±^) - \ Hv) - C(2) + 2 Li 2 (I±^) . 

In the considered case it is possible to obtain an explicit result in terms of hy- 
pergeometric functions and hence to perform the e expansion using directly the exact 
formulae. However, from our experience with the e expansion at d = 4 — 2e, we may 
say that this is a quite formidable and errorprone task. Thus, having done it once, it is 
much simpler to use the d- recurrence relations in order to get the result in the other 
dimension. In cases when the analytic result is not known, it might be useful to choose 
the most convenient d from the point of view of the e expansion and then to trans- 
form the result to the required space-time dimension. Especially this procedure maybe 
helpful if one encounters integrals with infrared as well as ultraviolet divergences. 

7. Conclusions 

We proposed a new type of recurrence relations for Feynman integrals. The addition 
of the space-time dimension to the set of recurrence parameters for Feynman integrals 
is quite natural and it extends the set of recurrence relations obtained from the inte- 
gration by parts method. We expect that within the framework of the extended set 
of recurrence relations it will be possible to formulate rather effective algorithms for 
the computation of Feynman diagrams. The problem with irreducible numerators in 
Feynman integrals finds a natural interpretation in the frame of generalized recurrence 
relations. We demonstrated that a direct solution of the new recurrence relations is 
possible. Finding solutions of such relations would considerably simplify if it would be 
possible to develop an effective algorithm for the asymptotic expansion at \d\ — > oo. 
At the same time the expansion at \d\ — > oo could be used as an approximation for the 
integral at finite d. This kind of expansion is well known in solid state and statistical 
physics. Basically, in the present publication we formulated general ideas about new 
techniques which we intend to apply later to the calculation of some particular classes 
of Feynman integrals. 
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